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Introduction 


•  Filtering  schemes  offer  a  less  dissipative  alternative  to  the 
standard  artificial  dissipation  operators  when  applied  to  high- 
order  spatial/temporal  schemes 

•  Ljmjtjng  Fact:  Filters  impart  the  same  amount  of  dissipation 
each  time  they  are  applied  so  if  they  are  applied  too  frequently, 
they  are  overly  dissipative 

•  Limiting,  Fact:  Stiff  systems  require  a  preconditioned  dual-time 
framework  to  be  solved  efficiently 

•  Limiting  Fact:  Filtering  cannot  be  applied  only  at  the  physical¬ 
time  level  and  convergence  guaranteed 

•  Limiting  Fact:  Filtering  cannot  be  straightforwardly  applied  at 
the  pseudo-time  level 

Objective:  T o  recast  common  filtering 
operators  as  equivalent  artificial-dissipation 

schemes 
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\w  Governing  Equations 


Navier-Stokes  Equations 

9Vi 


dCj  dFj  _ _ 

dt  dxi  dx: 


+  H 


Q  =  [p  pui  pe0]T 


Fj  —  'ypv.,  pu^u j  “t-  p&ij 

Quasi-linear  Form: 

9Q  5Q  =  avi 

dt  dx,  dx< 


T 

Uipho\  where  ho  =  eo  +  - 

P 


A  =  =  MAM"1 


A  =  diag  {ui  +  c,ui,ui  —  c} 


Primitive  Variables  Form: 

<9F,  av< 


P  <9QP 


dt  dxi  Ox, 


+  H 


dQp  dQp 

~e  dt  + - dx. 


dVj 

dxi 


+  H 


where  Te  =  and  Qp  =  [p  Ui  T]t 
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%#  Governing  Equations 


Residual  Form: 

+  Rs  (Q)  =  ^  +  Te  1Rs  (Qp)  =  0  where  Rs 

Temporal  Discretizations: 

-  SSPRK(3,  3): 

Q(1)  =  Q”  -  AtR,  (Qn) 
q(2)  =  5q»  +  I (q<D  _  AtRs  (qW)  ) 

Qn+1  =  ^Q"  + 1  (q(2)  -  AiRs  (q<2>)  ) 

-  SSPRK(4,  3):  Q(1)  =  Q 

Q(2)  =  Q(1)  -  ^AtR, 

q(3)  =  2Q„+1^Q(2)_  lAtR^Q(2)^ 

Q"+1  =  Q<3>  -  iA«Rs  (q<3)) 


dFj 

dxi 


dV, 

dxi 


-H 


Q 


Qp1*  =  Q"  -  AtRs,p  (Q") 

Qp2>  =  jQ?  + j(Qp* -^Rs,p(Qp1))) 

”+1  =  5Q;‘  +  5(Ql.2)-AiR5,p(Q‘2))) 

where  R„,p  =  rj1>wRs  (<#) 


'  2^fR«  (Q") 
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Spatial  Discretizations 


•  Discretely-conservative  High-order  Upwind  using 
Roe’s  Approximate  Riemann  Solver: 

^  =  F*"j+V2~  where  FiJ±1/2  =  F LR  =  f  (QL,  QR) 

oxi  A  Xi 

For  instance,  the  standard  second-order  central  difference  can  be  found  using: 
F+1/2  =  F(Q+0+F(Qo)  and  F_i/2  =  F(Qo)±F(Q-,) 

•  High-order  Interface  States: 


VL,+ 1/2 
VR,  +  1/2 
VL,  +  '/2 
VL,  + 1/2 
VL,+ 1/2 


-v-i  +  5n0  +  2v+i 

6 

-u_|_2  +  5n_)_i  +  2vq 
6 

2v-2  ~  13u_i  +  47uo  +  27n+i  —  3v+2 

60 

— 3n_3  +  25n_2  —  101u_i  +  319no  +  214v+i  —  38u+2  +  4u+3 

420 

4u_4  —  41n_3  +  199n_2  —  641u_i  +  1879uo  +  1375u+i  —  305u+2  +  55u+3  —  5u+4 

2520 
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Limiter:  Monotonicity-preserving 

Method 


Limits  the  left  and  right  interface  quantities  such  that 
they  are  within  an  interval  that  is  guaranteed  to 
preserve  monotonicity  with  an  appropriate  CFL  number 


Uses  a  five-point  stencil  to  distinguish  between  local 
extrema  and  discontinuities 


J  I  I  L 


J  L 


O 


J  I  I  l 


Further  details  are  in  the  paper 


G 


1 
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Limited  Quantities 


•  Either  the  conserved,  the  primitive,  or  the  characteristic 
variables  can  be  limited  by  the  MP  scheme 

*  Local  characteristic  limiting  can  be  carried  out  as  follows: 

For  Ql,+i/2!  the  stencil  consists  of  the  set  {Q_2,  Q-i,  Qo,  Q+i,  Q+2}- 

Wfc  =Mo1Qfc  f°r  *  =  -2,2 

Wl>+i/2  =  Mo  1Ql,+i/2 

or 


Wfe=N^Qp,fc  for  k  =  —2, 2 
WLj  +  1/2  =  N0  1Qp,L!  +  l/2 


where  N  1  =  M  1Te.  It  is  important  to  note  that  this  necessarily  implies  that 
the  same  Te  matrix  (calculated  at  the  midpoint  of  the  set)  resident  in  N"  _1  is 
used  for  all  five  points  in  the  limiting  stencil. 
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Sod’s  Shock  Tube  Problem 


•  Domain  20  m  long 

*  100  points 

•  Discontinuity 
initially  at  x  =  10  m 

*  Run  for  50  time 
steps 


(a)  Characteristics  from  Conserved 


(b)  Characteristics  from  Primitive 


(c)  Conserved  (d)  Primitive 
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Lax’s  Problem 


Domain  20  m  long 

100  points 

Discontinuity 
initially  at  x  =  10  m 

Run  for  160  time 
steps 


PL 

UL 

PL 


At  =  0.0002s 
cr  =  0.403 
if  —  0.032s 


0.890,  pR  =  1.000 
69.8,  Ur  =  0.0 

70560,  pR  =  11420 


(a)  Characteristics  from  Conserved 


(b)  Characteristics  from  Primitive 


(c)  Conserved  (d)  Primitive 
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1 D  Contact  Discontinuity  Problem 


Domain  20  m  long 

100  points 

Discontinuities 
initially  at  x  =  5  m 
and  x  =  15  m 

Run  for  16,000  time 
steps 


1100 

1000 

900 


600 


300 

200 


1 

! - 1 - 

— i - 1 - 1 - 1 - 1 — 

- 1 - 1 - 

- 

G 

3 

C 

i  i  j  i 

3 

. 

- 

3  j  C 

- 

3  G 

_ i _ i _ 

10 

X 


12  14  16  18  20 


(a)  Characteristics  from  Conserved 


At  =  6.25  x  10~5s 


(T  =  0.370 
tf  =  1.00s 
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(c)  Conserved 
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(b)  Characteristics  from  Primitive 


(d)  Primitive 
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Temperature  Temperature 


ID  Contact  Discontinuity  Problem 


(a)  Characteristics  from  Conserved 


(a)  a  =  4 


(b)  a  =  7 


(b)  Characteristics  from  Primitive  (c)  a  =  1 1  (d)  a  =  1 9 
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2D  Drop  Problem 


•  Domain  1  mm  square 

•  100  points  square 

•  Cold  Fluid  inside  the  drop 

•  Run  for  25,000  time  steps 

At  =  4.0  x  10-9s 
a  ~  0.4 

tf  =  1.0  x  10-4s 

PL  =  15,  pR  =  1258 
Ul  =  10.0,  Ur  =  10.0 
Pl  =  1  x  107  pr  =  lx  107 


(a)  Density,  Exact  (b)  Temperature,  Exact 


(c)  Density,  Characteristics  from  Conserved  (d)  Temperature,  Characteristics  from  Conserved 


Q  =  Q rF  +  Ql  (l  —  F) 


1  —  tanh 


r  -  r_d 

^ d 


r d  =  0.2mm  and  =  0.015mm 


(e)  Density,  Characteristics  from  Primitive  (f)  Temperature,  Characteristics  from  Primitive 
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Temperature  Temperature 


2D  Drop  Problem,  Centerline  T 


(a)  Characteristics  from  Conserved 


(a)  a  =  19 


(a)  Characteristics  from  Conserved,  2x  Resolution 


(b)  Characteristics  from  Primitive 


(b)  a  =  19,  Ninth-order 


(b)  Characteristics  from  Primitive,  2  x  Resolution 
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PERIODIC 


2D  Rotating  Detonation  Wave 


•  Unwrapped  cylinder 

•  Follows  the  work  of 
Schwer  (reference  in 
paper) 


(a)  Limit  All  Primitives 


T  [K] 


13075. 
2369. 
1662. 


956.2 

^250.0 


•  Stoichiometric  H2/Air 

*  P plenum  =  10  ATM 
T plenum  =  300K 


T  [K] 


-956.2 

^250.0 


OUTLET 


(b)  Limit  Non-monotone  Primitives 


T[K1 
1-3075. 


2369. 

1662. 


956.2 
*  250.0 


(c)  Limit  Non-Monotone  Characteristics  From  Primitives 
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Conclusions 


•  Applying  existing  high-order  monotonicity¬ 
preserving  methods  to  the  primitive  variables  is 
not  straightforward 

•  The  MP  scheme  should  be  applied  to  the 
characteristic  variables  calculated  from  the 
conserved  variables  for  the  best  results 

•  These  results  will  offer  unique  results  for  more- 
complicated  equations  of  state  (EoS)  where 
temperature  and  pressure  cannot  be  calculated 
from  density  and  internal  energy  directly 
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Future  Work 


•  Use  a  detection  strategy  to  limit  the  frequency  with 
which  the  limiter  and,  with  it,  the  conversion  from 
conserved  to  primitive  via  the  EoS  is  applied 

-  Interpolate  primitive  variables 

-  Detect  where  limiting  might  be  needed 

-  Only  where  needed,  calculate  characteristic  variables 
from  conserved  variables 

-  Note  specifically  where  the  MP  limiter  is  active 

-  Convert  back  to  primitive  variables  only  where  the  limiter 
is  active 

-  Use  the  first-order  accurate  value  for  the  initial  guess  of 
temperature  and  pressure  in  the  EoS  routine 

•  In  this  way,  efficiency  and  accuracy  can  be  optimized 
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Questions??? 
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